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摘要 :为 了 确定 月 球 与 即将 发 射 到 地 月 系 La 点 的 中 继 星 之 间 的 位 置 关 系 进 而 估算 中 继 星 激光 测 


距 的 成 功率 ， 按 照 嘿 轨道 周期 为 14 天 左右 的 要 


求 对 中 继 星 所 在 的 晕 轨 道 进行 了 计算 ， 并 建立 了 一 人 


综合 考虑 望远镜 抖动 、 大 气 抖动 和 预报 轨道 横向 偏离 的 模型 。 从 数值 上 给 出 了 一 条 轨道 周期 为 


14.78 R, XA (地 月 连 线 方向 ) 的 振幅 为 12493km， 了 方向 为 34596km，Z 方 向 (垂直 于 地 月 轨道 


平面 方向 ) 为 11916km 的 周期 轨道 。 由 于 坚 轨 道 的 最 小 振幅 远大 于 月 球 侦 挡 的 临界 振幅 4000km， 因 


此 月 球 对 中 继 星 不 存在 侦 挡 问题 。 基 于 上 面 建立 的 测 距 成 功率 模型 ， 根 据 昆 明 站 (国际 编号 : 7820) 


的 激光 测 距 系 统 对 运行 在 该 尝 轨 道上 的 中 继 星 进行 测 距 成 功率 分 析 ， 结 果 表明 : 测 距 成 功率 随 着 中 


继 星 横 向 轨道 标准 差 的 增 大 呈现 快速 降低 的 趋势 。 对 于 中 继 星 到 测 站 的 平均 距离 而 言 ， 当 中 继 星 没 


有 横向 偏离 时 ， 探 测 器 产生 的 光电 子 数 为 0.151， 成 功率 为 14.07%; 横向 偏离 2km 时 ， 光 电子 数 降 


为 0.035， 成 功率 降 为 3.46%。 对 比 最 近 距 离 与 最 远 距离 的 情况 ， 无 横向 偏 移 的 情况 下 ， 探 测 器 产生 
的 光电 子 数 从 0.174 J 0.139， 成 功率 从 16.01% 降 为 13.02%。 根 据 上 面 所 有 的 计算 结果 ， 为 后 续 


云南 天 文 台 1.2m 望远镜 实现 中 继 星 激光 测 距 提 


供 参 考 。 
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引力 波 的 探测 不 仅 是 检验 广义 相对 论 根基 的 重要 途径 ， 同 时 也 为 研究 宇宙 提供 了 一 种 全 新 的 手 


段 。 继 地 基 的 LIGO 和 Virgo 同时 直接 探测 到 
引力 波 探测 计划 也 在 紧锣密鼓 地 进行 中 ， 其 中 ， 
近 圆 轨道 上 通过 观测 卫星 之 间 的 距离 变化 探测 引 


105km， 轨 道 平面 面向 由 两 颗 快 速 绕 转 的 白矮星 


“天 琴 计 划 ” 的 第 1 步 是 实现 激光 测 月 和 


文章 编号 :1672-7673(2019) 


由 双 中 子 星 合并 产生 的 引力 波 和 伽 马 射线 目 ， 天 基 的 

“天 琴 计 划 ” 是 将 3 个 星 载 激光 干涉 仪 发 射 到 绕 地 
| 力 波 外 。 这 3 颗 卫 星 组 成 等 边 三 角形 构 型 ， 辟 长 约 为 
组 成 的 引力 波 辐射 源 RXJ0806.3+15276]。 
深 空 卫星 激光 测 距 。 为 了 满足 即将 发 射 的 “ 娣 娥 4 


号 ”月 球 背面 登陆 器 与 地 球 的 实时 通讯 ， 需 要 地 


E 地 月 系 L2 点 附近 的 坚 轨 道上 放置 一 个 中 继 星 实 现 通 


讯 中 转 外 。 文 [5] 针 对 地 月 工 ,点 的 中 继 星 任 务 分 析 了 发 射 窗口 、 地 月 转移 时 间 、 近 月 点 高 度 、 近 月 点 
倾角 、 轨 道 振幅 等 因素 对 地 月 Ls 转 移 轨 道 和 使 命 轨 道 特性 的 影响 ， 文 [6] 针 对 中 继 星 转移 轨道 的 转 


移 时 间 、 近 月 点 高 度 和 尝 轨 道 振 幅 等 约束 条 件 ， 


究 了 地 月 系 LL 点 纯 反射 式 激光 测 距 技术 和 任务 
给 出 了 预期 系统 能 接收 到 的 单 脉 冲 回 波光 子 数 。 


OE THREF L 点 转移 轨道 设计 方法 。 文 [7] 研 
设计 ， 并 结合 云南 天 文 台 1.2m 望 眼镜 激光 测 距 系统 
由 于 中 继 星 的 星 等 (上 暗 于 17 等， 计算 过 程 见 附 录 ) 


高 于 1.2m 望远镜 光电 探测 的 极限 星 等 ， 所 以 有 必要 计算 出 一 条 假想 的 晕 轨 道 并 对 预报 精度 如 何 影 
响 单 脉冲 回 波 光子 数 做 出 分 析 ， 进 而 研究 评估 中 继 星 和 月 球 处 于 不 同位 置 关系 时 噪声 对 激光 测 距 成 
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功率 的 影响 。 本 文 根 据 轨 道 对 称 性 理论 和 微分 改正 法 馈 计 算出 一 条 围绕 地 月 系 L 点 运动 的 坚 轨道 ， 
在 此 基础 上 估算 了 能 量 服从 高 斯 分 布 的 光斑 在 考虑 大 气 抖动 、 望 远 镜 抖 动 、 中 继 星 预 报 精度 的 情况 
下 ， 单 脉冲 的 回 波光 子 数 和 测 距 成 功率 随 中 继 星 预报 精度 的 变化 。 
1 围绕 L 点 的 党 轨道 计算 方法 

中 继 星 在 地 月 系 中 的 圆 形 限 制 性 三 体 问 题 ( Circular Restricted Three Body Problem, CRTBP) 指 
的 是 地 球 和 月 球 绕 地 月 系 质心 作 圆 周 运动 ， 中 继 星 在 地 球 和 月 球 共同 引力 作用 下 的 动力 学 问题 。 这 
里 仅 考虑 中 继 星 围绕 L 点 的 周期 运动 ， 并 最 终 给 出 一 条 可 行 的 晕 轨 道 。 

定义 地 月 距离 4 为 基本 的 长 度 单位 ， 地 月 周期 P 为 基本 的 时 间 单 位 ， 那 么 无 量 纲 的 距离 单位 
Lx*、 时 间 单 位 tk、 速 度 单 位 vy* 与 正常 使 用 的 距离 单位 元 、 时 间 单 位 t+、 速度 单位 vy 之 间 的 转换 关 
系 如 下 : 


L*=Lid , 
t*=2r*t/P , 
v*= Lt v*P/2*n*d (1) 
其 中 : 
(2) 
其 中 ，G 为 引力 常数 ，M 为 地 球 月 球 质 量 之 和 。 圆 形 限 制 性 三 体 问题 在 绕 质 心 旋转 坐标 系 下 的 无 量 
岗 运动 方程 站 为 


x Ver 

y vy 

之 Uz 

Var 2v, + ge 
Da, au 

Uy 2Ux 4 By 

Uv au 

‘2 Oz 


(3) 


Uyd =(1/2)*(02+ 2) + (p/n + ur, (4) 
r=[e+ py ty +z)” 
m=) ty te] , (5) 
KH=m/(mtm) . (6) 
(ORF, m 和 m 分 别 为 地 球 和 月 球 的 质量 。 圆 形 限 制 性 三 体 问题 存在 一 个 Jacobian 积分 ， 即 
2UGep2) Ht) =C (7) 
其 中 C 为 一 个 常数 。 如 果 令 中 继 星 的 速度 为 零 ， 那 么 会 得 到 零 速度 曲面 2UG,y,z) = C， 更 进 一 
步 ， 令 z=0， 可 以 得 到 零 速度 曲面 在 系 了 平面 的 投影 ， 即 零 速度 曲线 2UCcy) = C。 
在 加 时 刻 ， 取 一 个 位 于 Lo 点 附近 的 初始 状态 矢量 (xo,0,zo0,vyo0)， 那 么 经 过 半 个 周期 tc 后 中 继 
星 的 状态 矢量 可 以 表示 为 
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x(to + T) $1(%0, 0, 20, 0, vo, 9, T) 
y(to +T) $2(20, 0, 20,9, Vy), 9, T) 
z(to +7) 3 To “0,2,0, fae T) 


£o, 0, 20, 0, vyo, 0,7) 


’ yo 0 


vz(to + T) Qa 
vy(to + T) s(x, 0, 20, 0, w0, T) 
vz(to 十 T) 6(Xo, 0, 20, 0, Vy 9, T) 


( 
( 
_ | ¢a( 
( 
( 
( 


(8) 
EA, :(%0,0,20,0,vy0,0,t) 为 一 个 流 。 从 初始 位 置 经 过 一 个 周期 2rz 后 ， 位 置 偏 量 定义 为 


Ar =((x(to + 2t)-X0) + (y(to + 2t)-yo) + (z(to + 21)-2)?)!” (9) 
根据 对 称 性 原理 中 ， 经 过 半 个 周期 + 后 中 继 星 的 状态 矢量 满足 


po (To, 0, Zo, 0, vyo, 0, 7) 
d(x, 0, 29,0, Vya; 9, T) =0 
$6(Xo, 0, 20, 0, Vy, 0, T) 
(10) 
固定 xo M C10) 式 含 有 3 个 未 知 数 ， 所 以 采用 微分 改正 法 求 zo, vyor: 


J Ezz p 
Then Lra |28 k Bho [Pl ae 
w nt aD 
其 中 ,6gxy6zo= 5B(2,3) 
Ob »/vyo=D(2,5) 
06 »/Ot= v, 
Op 4/Oz0=B(4,3) 
Ob s/Ovyo= (4,5) 
Op 4/Ot= 2v, +OU/Ox 
Ob /Oz0= (6,3) 
Ab</Ovyo= B(6,5) 
Ob</Or=OUI0z (12) 
H i en (关于 状态 转移 的 定义 和 推导 见 附录 ) 。 当 nK OEF, [zo,vy0,t]'0 toH} 
刻 的 初始 值 ， 经 过 一 次 微分 改正 后 ， 得 到 一 个 新 的 初始 值 [zovyor] i » ERRAR (Cn 足够 大 ) 会 得 到 


满足 精度 要 求 的 初始 值 [zovyo,t]%。 对 于 Li 点 用 上 述 方法 可 以 得 到 一 条 三 维 晕 轨 道 ， 而 对 于 LA, 
用 这 种 方法 只 能 得 到 平面 上 的 Lyapunov $k, AA a MEMEO. AT RSS BHI, 
用 前 面 得 到 的 平面 Lyapunov 轨道 作为 初始 条 件 。 固 定 zo， 使 xo,vyo，7 为 变量 ， 那 么 《10) REF 3 
个 未 知 数 ， 因 此 采用 微分 改正 法 来 求 xo vyo t: 
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， 其 中 

cqpvoxm= (2, 1) 
Op/0vo= O (2, 5) 
b2/d T= vy 
Op/0%= oO (4, 1) 
aqpyavo= O (4, 5) 
0/0 T= 2v, +0U/ðx 
OM/Ox0= (6, 1) 
OM/Ovyp= 0 (6, 5) 
Om /0 T=0U/0z 


以 Li 为 坐标 原点 ， 选 取 初始 位 置 为 (74524.0,0,-2000.0)km， 初 始 速 
半 周 期 为 6.30 天。 通过 〈1) 式 和 “(2 式 ) 进行 无 量 纲 化 后 ,代入 运动 方程 ， 并 积 
个 平面 Lyapunov 轨道 的 初始 条 件 

该 轨道 周期 为 
EJ 1.57801632024x10"', Jacobian 积分 常数 为 3.15056044173， 状 态 转 


转移 矩阵， 求解 微分 改正 方 


14.8485511785 天 ， 位 置 偏 量 


Op» 


Oxo 
O64 


Odo 


Or 


Obs 


“BO Oro 


T ape 


Oxo 


(14) 


Or 


Ope 


OT 


程 (11)， 最 后 得 到 满足 
(1.1817143086500759,0,3.5933035277813563 10-22,0,-0.16170712205794957,0)。 


n|to+rT 
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(13) 


KEEN (0,-2763.0,0)km/day， 和 初始 


运动 方程 和 状态 


移 和 矩阵 的 模 为 1.00000000018 (接近 于 1) ， 说 明 该 轨道 满足 周期 轨道 条 件 。 将 平面 Lyapunov 轨道 


的 初始 条 件 做 一 个 小 的 改动 ， 
运动 方程 和 状态 转移 矩阵 ， 
曼 轨 道 的 2 继续 做 一 个 小 的 
求解 微分 改正 方程 (13) ， 


例如 将 z 改动 到 1.5 x 107"， 其 他 不 变 。 


求解 微分 改正 方程 
改动 ， 其 他 不 变 ， 
不 断 重复 。 最 后 得 到 


13) ， 


将 其 代入 运动 方程 ， 并 积分 


二 


得 到 一 个 稍微 远离 平面 的 尝 轨 道 ， 然 后 将 该 


再 次 代入 运动 方程 ， 


并 积分 运动 方程 和 状态 转移 矩阵 ， 


一 条 轨道 周 为 14.7843020586 天 ， 初 始 条 件 为 


(1.179549767505286,0,0.03662109375,0,-0.16319295932416145,0) 的 曙 轨 道 〈 见 图 


的 位 置 偏 量 为 1.5556312755 


9x 107, Jacobian 积分 
0.999999999964 (接近 于 1) ， 说 明 满 足 周期 轨道 条 件 。 


1 和 图 2) 。 该 轨道 
常数 为 3.14635368089， 状 态 转 移 矩 阵 的 模 为 
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1 地 月 系 工 ,点 的 尝 轨道 三 维 轨迹 图 及 其 在 3 个 平面 的 投影 。 图 中 的 尺度 单位 为 地 月 距离 ， 即 384402km， 月 
球 大 小 和 党 轨道 尺度 为 等 比例 画 出 ， 原 点 位 于 地 月 系 质心 。 


Figl. Three-dimensional trajectory of Earth-Moon halo orbits around L2 and its projection on three 


orthogonal planes. The unit of length scale is in distance between Earth and Moon, namely 384402km; the 
size of Moon and halo orbit are plotted in equal proportion, and the origin is located at the mass center of 


Earth-Moon system. 
IKPaR AE 2 UE YZ PIMP ay, UA A RD AEA 1.24, mH PRAY 
角 半 径 约 为 0.26， 因 此 满月 时 月 球 的 光 背 景 可 能 对 中 继 星 激 光 测 距 产生 一 定 影响 。 从 文 [6] 的 计算 
可 知 ， 当 晕 轨 道 的 振幅 不 小 于 4000km 时 即 可 实现 月 球 无 遮挡 ， 本 文中 时 轨道 的 最 小 振幅 为 
23832km， 因 此 月 球 对 该 尝 轨 道 没 有 遮挡 。 在 确认 满月 对 中 继 星 激光 测 距 的 影响 后 ， 下 面 对 中 继 星 
的 测 距 成 功率 展开 建 模 与 分 析 。 


0.3 Halo orbit projection and Jacobian Constant in X-Y plane 


0.2 


0.1 


> 0.0 L1 Moon 


—0.1 


—0.2 


“0.7 0.8 0.9 1.0 1.1 1.2 
x 


图 2 地 月 系 L RABEL AE X- 了 平面 的 投影 和 零 速 度 曲 线 


Fig2. Projection of the Earth-Moon halo orbit around L2 and its zero velocity curve in X-Y Plane. 


2 测 距 成 功率 模型 
测 距 成 功率 指 的 是 单 光 子 探测 器 能 够 探测 到 至 少 一 个 光电 子 的 概率 。 由 于 探测 器 探测 光电 子 的 
概率 服从 泊 松 分 布 "3， 即 
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P(k:2) = e-A* (Ak) , (15) 
测 距 成 功率 可 以 表示 为 


= 


¢= 1-P(0;N) = 1-e (16) 

其 中 , N 为 探测 器 产生 的 平均 光电 子 数 。 光 电子 数 越 多 信号 越 强 ， 反 之 光电 子 数 越 少 信号 越 弱 。 计 算 

测 距 成 功率 的 关键 是 计算 探测 器 产生 的 平均 光电 子 数 ， 下 文 主要 围绕 如 何 计算 探测 器 产生 的 平均 光 

电子 数 展开 讨论 。 
在 不 考虑 望远镜 抖动 、 大 气 抖 劲 和 中 继 星 偏离 预报 轨道 的 理想 情况 下 ， 望 远 镜 发 射 的 激光 光斑 

中 心理 论 上 应 该 对 准 中 继 星 ， 此 时 激光 光斑 的 能 量 密度 服从 高 斯 分 布 (Gaussian Distribution) , El 

, (17) 
其 中 , Eo 为 光斑 的 总 能 量 ，p 为 光斑 上 任意 一 点 到 光斑 中 心 的 距离 ，pe 为 光斑 的 特征 半径 并 满足 如 下 
关系 式 : 


S 


一 


. (18) 
取 光 斑 的 半径 为 1.5 倍 特征 半径 时 ， 光 斑 内 的 能 量 占 总 能 量 的 98.89%。 光 斑 半 径 r 与 激光 发 散 角 
0, Hee SPAS ER RS R 以 及 望远镜 有 效 口径 DD 之 间 的 关系 为 

r=1/2*(RO+D) (19) 
HARRES KB) 以 及 中 继 星 的 轨道 偏离 时 ， 光 斑 的 平均 能 量 密 度 分 布 可 以 表 
RA 


(20) 
其 中 ,望远镜 抖动 g(x,y) ~M0,0D; KARS g(x,y) NOAD: 预报 轨道 偏离 gay) NO; J 
为 三 阶 单位 矩阵 ， 积 分 区 域 为 整个 光斑 覆盖 面 RR; * 表示 卷 积 运算 叫 。 对 (20) 式 进行 推导 并 化 简 得 


(21) 

大 气 抖动 引起 光斑 中 心 漂移 的 均 方差 os 表示 为 5 
(22) 

其 中 ;nm 为 大 气相 干 长 度 ; k= 20/1, A 为 激光 波长 。 探 测 器 探测 到 的 光电 子 数 ON 可 以 表示 为 "5 
(23) 


EP, 为 望远镜 发 射 系统 的 光学 效率 ， 芳 为 大 气 透 过 率 ; LAER RO, SA PARRA 
反射 面积 ;peeeion 为 反射 器 的 反射 率 ; assio 为 中 继 星 反 射 器 的 反射 立体 角 ， 六 为 望远镜 接收 系统 
的 光学 效率 ;ma 为 探测 器 的 量子 效率 ; h AELE e 为 真空 中 的 光速 。 中 继 星 反 射 器 的 反射 立 
体 角 Quettction 与 反射 发 散 角 Protection 的 关系 如 下 


(24) 
3 结果 与 分 析 
在 中 继 星 激光 测 距 过 程 中 ， 取 激光 脉冲 的 能 量 Eo 为 3000mJ， 测 站 到 中 继 星 的 平均 距离 R 取 测 
站 到 工 点 的 距离 ， 即 442548km， 望 远 镜 的 抖动 标准 差 of A 0.1”( 乘 以 R 换算 成 长 度 单位 ) ， 大 


Stephenson P C. Satellite laser ranging photon-budget calculations for a single satellite cornercube retroreflector: attitude 
control tolerance 
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气 的 相干 长 度 to 10cm, W 
射 器 的 反射 率 Preflection 为 0.0， 
TH 0.1， 望远镜 发 射 系统 的 
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光 发 射 的 发 散 角 9 为 2"”， 中 继 星 的 有 效 反射 面积 8 为 0.0227m’, fz 
反射 器 的 发 散 角 Qienection 为 之 ， 大 气 透 过 率 Ta H 0.6, 卷 云 透 过 率 
光学 效率 六 为 0.4， 接 收 系统 的 光学 效率 六 为 0.2， 探 测 器 的 量子 效 


率 四 为 0.6， 望 远 镜 的 有 效 口径 D 为 1.06m， 激 光 的 波长 4 为 532nm07。 通 过 前 面 计 算 的 又 轨道 


可 知 ， 测 站 到 中 继 星 的 最 近 距 


离 和 最 远 距 离 分 别 为 427287km 和 451889km ， 保 持 其 他 参数 不 变 ， 


预报 轨道 横向 标准 差 oH 0 ~5km 时 探测 器 产生 的 光电 子 数 和 测 距 成 功率 见 图 3。 
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图 3 GETZ GER) 和 测 距 成 功率 〈 实 线 ) 随 中 继 星 的 横向 预报 轨道 标准 差 的 变化 图 ， 从 左 到 右 分 别 对 应 测 


到 中 继 星 最 近 距 离 、 平 均 距离 、 最 远 距 离 的 情况 


Fig3. Variations of the photoelectron numbers(dotted line) and the laser ranging success probability(solid 


line) with the transverse standard deviation of the predicted relay satellite orbit; from left to right are 


situations for the nearest, average, furthermost distance from station to the relay satellite. 


从 图 3 中 可 看 出 光电 子 数 和 测 距 成 功率 随 着 中 继 星 轨道 横向 标准 差 的 增 大 而 迅速 减 小 ， 尤 其 是 
横向 标准 差 位 于 1 ~2km 时 下 降 最 快 ， 超 过 2km 后 下 降 变 组 。 对 于 平均 距离 而 言 ， 当 中 继 星 没 有 偏 


离 预报 时 ， 探 测 器 产生 的 区 


EF% 0.1$1， 成 功率 为 14.07%， 偏 离 2km 时 的 光电 子 数 降 为 


0.035， 成 功率 降 为 3.46%， 此 时 光电 子 数 和 成 功率 分 别 衰减 了 76.8% 和 75.4%; 对 于 最 近 距 离 而 


言 ， 当 中 继 星 没有 偏离 预报 时 


， 探 测 器 产生 的 光电 子 数 为 0.174， 成 功率 为 16.01%， 偏 离 2km 时 


的 光电 子 数 降 为 0.038， 成 功率 降 为 3.77%， 此 时 光电 子 数 和 成 功率 分 别 衰减 了 78.2% 和 76.4%; 


对 于 最 远 距离 而 言 ， 当 中 继 


星 没有 偏离 预报 时 ， 探 测 器 产生 的 光电 子 数 为 0.139， 成 功率 为 


13.02%, WA 2km 时 的 光电 子 数 降 为 0.033， 成 功率 降 为 3.29%， 此 时 光电 子 数 和 成 功率 分 别 衰减 


了 76.2% 和 74.7%。 通 过 对 上 
0.139， 成 功率 从 16.01% 降 为 


上 最近 距 离 与 最 远 距 离 的 情况 ， 探 测 器 产生 的 光电 子 数 从 0.174 降 为 
13.02%. 


oA ”Vi 人 人 [人生 甘 日 工 | 上 | 
CI I law IV | F HH T lJ 


4 结语 
本 文 针对 即将 发 射 到 地 月 系 Lo 点 的 中 继 星 进行 了 坚 轨 道 计 算 ， 根 据 轨 道 的 对 称 性 理论 和 微分 

改正 法 ， 从 数值 上 给 出 了 一 条 轨道 周期 为 14.78 R, XTRA (地 月 连 线 方向 ) 的 振幅 为 

12493km， 了 方向 为 345$96km，Z 方 向 (垂直 于 地 月 轨道 平面 方向 ) 为 11916km 的 轨道 。 运 行 在 该 轨 

道上 的 中 继 星 与 月 球 的 角 距离 约 为 1.24， 而 月 球 的 角 半 径 约 为 0.26， 因 此 月 球 的 光 背 景 对 中 继 星 

测 距 可 能 会 产生 一 定 的 影响 。 另 一 方面 ， 由 于 轨道 的 最 小 振幅 远大 于 月 球 遮挡 的 临界 振幅 
4000km， 所 以 月 球 对 中 继 星 不 存在 遮挡 问题 。 中 继 星 在 V 波段 的 视 星 等 暗 于 17 等 ， 已 经 超出 了 
1.2m 望远镜 的 极限 星 等 ， 况 且 满 月 时 望远镜 受 月 光 的 影响 ， 更 加 增加 中 继 星 可 见 的 难度 。 随 后 在 综 
合 考虑 大 气 抖动 、 望 远 镜 抖 动 和 中 继 星 横向 偏离 预报 轨道 的 情况 下 ， 建 立 了 一 个 估算 回 波光 电子 数 
和 测 距 成 功率 的 模型 。 根 据 云 台 的 激光 测 距 系统 对 运行 在 该 尝 轨 道上 的 中 继 星 进 行 回 波光 电子 数 和 
测 距 成 功率 分 析 ， 结 果 表明 光电 子 数 和 测 距 成 功率 随 着 中 继 星 轨道 横向 标准 差 的 增 大 而 迅速 减 小 ， 

™ 尤其 是 横向 标准 差 位 于 1 ~2km 时 下 降 最 快 ， 超 过 2km 后 下 降 变 缓 。 对 于 中 继 星 到 测 站 的 平均 距离 

i 而 言 ， 当 中 继 星 没 有 横向 偏离 时 ， 探 测 器 产生 的 光电 子 数 为 0.151， 成 功率 为 17.07%; 横向 偏离 

2km 时 ， 光电 子 数 降 为 0.035， 成 功率 降 为 3.46%。 对 比 最 近 距 离 与 最 远 距 离 的 情况 ， 无 横向 偏 移 

的 情况 下 ， 探 测 器 产生 的 光电 子 数 从 0.174 降 为 0.139， 成 功率 从 16.01% BEN 13.02%。 尺 管 最 远 

距离 与 最 近 距 离 相差 24602km， 但 探测 器 产生 的 光电 子 数 没有 明显 的 改变 ， 这 是 由 于 轨道 的 振幅 相 

对 于 整个 地 月 距离 而 言 较 小 的 缘故 。 

附录 

1 状态 转移 矩阵 的 定义 和 推导 

给 定 一 个 动力 学 系统 dX(b/dt=FCX() 和 to 时 刻 的 一 个 初始 条 件 X(to) 以 及 tb 时 刻 的 一 个 微小 
扰动 5X( tb)， 随 着 时 间 的 演化 ，t 时 刻 该 动力 学 系统 的 扰动 变 为 XA ME 4) o REFET E RE 
P(t, to) 描述 了 tb 时 刻 的 微小 扰动 SX(t) 与 t 时 刻 的 微小 扰动 sx) 之 间 的 关系 ， 即 


O D(t, to) =OX(t)/OX( to)» 
状态 转移 矩阵 实际 上 是 微分 方程 d@(t,to)/dt = 6X(t)/0X(t)*@(t, to) 的 解 ， 该 微分 方程 的 推导 如 下 : 
因为 


dX(t)/dt= F[X(t)], 
所 以 
d(X(t) + 8X(t)) dt= F[X(t) + 8X(t)], 


一 阶 展 开 后 ， 得 到 


dX(t) dt+d[8X(t)] dt= FLX(t)] +OF/OX*8X(t) + O[SX(D]， 
略 去 高 阶 项 并 化 简 得 
d[5X(t)]/dt=OF/OX*8X(t), 
由 于 
5X(D = O(t, to)*SX(to)， 
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将 其 代入 上 面 的 方程 得 到 
d@(t, to /dt*8X( to) =OF/AX*O 
最 后 化 简 得 到 
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(t, to )SX( to); 


A(t, to )/dt=OF/OX* W(t, to), 


对 该 微分 方程 积分 即 可 得 到 t RASA EE @B(t, tb)。 由 定义 可 知 to 时 刻 的 状态 转移 矩阵 


@(to,to)=I， 其 中 I 为 单位 矩阵 。 对 于 周期 轨道 有 O( to+27,to)=I， 其 中 1 为 半 周 期 。 


2 HRE V 波段 视 星 等 的 估算 
把 太阳 视 为 表面 温度 为 5778K 的 黑体 ， 那 么 单位 面积 站 
定律 


位 波长 的 太阳 辐射 功率 满足 普 朗 元 辐 冉 


图 4 微小 扰动 随时 间 的 演化 示意 图 


Fig4. Evolution of the perturbation with the time by a trivial perturbation on the initial value. 


其 中 ， 工 为 太阳 表面 温度 ;和 为 波长 ;bh 为 普 朗 克 常 数 ，c 为 光速 。 因 为 V 波段 的 光谱 范 
围 为 507nm ~595nm， 中 心 波 长 为 551nm， 因 此 对 MMT) 在 V 波段 积分 即 可 得 到 太阳 表 
面 单 位 面积 的 辐射 功率 Fs， 


其 中 ， MA 入 为 V 波段 的 下 限 波长 和 上 限 波 长 。 中 继 星 接收 到 的 太阳 辐射 功率 Fe 可 表 


示 为 


2 


其 中 ， Rs 为 太阳 的 半径 ，Rs 为 太阳 到 地 月 质心 的 距离 ，A 为 太阳 帆 的 展开 面积 ， 此 处 取 
值 为 10m。 由 中 继 星 太阳 帆 反 射 到 达 测 站 的 单位 面积 上 的 太阳 辐射 功率 F 为 


? 


KH, p 为 太阳 帆 的 反射 率 ， 一 般 取 值 为 16; R 为 中 继 星 到 测 站 的 平均 距离 ， 取 值 为 L 
点 到 测 站 的 距离 ， 即 442548km。 为 了 计算 中 继 星 在 V 波段 的 视 星 等 ， 需 要 一 个 基准 值 ， 
该 基准 值 可 取 视 星 等 为 0 的 辐射 源 。 测 站 单位 面积 接收 到 的 辐射 功率 Fu 可 表示 为 


更 


其 中 ， vi 和 vo 表示 V 波段 的 下 限 频率 和 上 限 频 率 ; Fo 为 单位 面积 单位 频率 的 功率 ， 取 
为 3640Jy。 最 后 根据 星 等 计算 公式 


m 


来 求 得 中 继 星 在 V 波段 的 视 星 等 ， 结 果 为 17.2. 
3 计算 学 轨 道 的 python 代码 
本 文 涉及 的 Python 代码 均 可 在 smellysheep.com 上 下 载 ， 主 要 包括 : 
。 计 算 地 月 系 5 个 平 动 点 的 位 置 ， 计 算 Jacobian 积分 常数 ， 积 分 运动 方程 和 状态 转 
移 矩 阵 ， 求 解 微分 改正 方程 ， 寻 求 L 点 的 平面 Lyapunov 轨道 和 三 维 晕 轨道 
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Estimates of the success probability for relay satellite laser ranging 
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Abstract: In order to identify the positional relationship between the Moon and the 
incoming relay satellite on a halo orbit around L+- of the Earth-Moon system, and then to 
evaluate the laser ranging success probability, in this paper we make a computation on 
the halo orbit with a request of about 14 days in its period, and establish a model given 
the tremble of the telescope, of the atmosphere and the transverse deviation from the 
predicted orbit. A numerical halo orbit is given with the period of 14.78 days, the 
amplitude of 24986km in X(along the line of Earth and Moon), 69192km in Y, 23832km in 
Z(orthogonal to the orbit plane of Earth and Moon). The minimum amplitude of the 
calculated halo orbit goes far beyond the Moon-sheltered critical amplitude which is 
generally 4000km, thus no shade from the Moon on the relay satellite is existed. Based 
on the model of laser ranging success probability established previously, we make an 
analysis on the success probability according the laser ranging system at KUNL 
station(International ID:7820), the results indicating that the laser ranging success 
probability rapidly drops down with the increase of the orbit transverse standard 
deviation. For the average distance from the station to the relay satellite, the detector 
produces about 0.151 photoelectrons within a single pulse and the corresponding 
success probability is 14.07% when there is no deviation from the predicted orbit; the 
number of photoelectrons drops down to 0.035 and the corresponding success 
probability down to 3.46% when the relay satellite transversely deviates 2km from the 
predicted orbit. Comparing between the nearest distance and the furthermost distance, 
with no deviation, the number of photoelectrons declines to 0.139 from 0.174 and the 
corresponding success probability down to 13.02% from 16.01%. All of the above results 
provide a reference for the realization of the relay satellite laser ranging with 1.2m 


telescope at Yunnan Observatory in the following days. 
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